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METHOD AND APPARATUS FOR CALCULATING THE REMAINDER OF A 

MODULO DIVISION 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates to a method and apparatus for calculating the 
remainder of a modulo division, and more particularly, to a non-iterative technique for 
5 calculating the remainder of a modulo division. 

DESCRIPTION OF THE RELATED ART 

Most computer systems require the value modulo n of an integer m, generally 
written m mod n. M mod n operations are utilized in a number of computer processes, for 
i/f example, address generating, communication channel load balancing, computer graphics, 

M= 10 telephone switching and telephone transmission, packet switching and transmission, and 
[ digital message encoding/decoding. Therefore, a routine that accepts integers m and n and 

^ produces m mod n is a useful computer tool. 

PJ The traditional method of calculating the remainder of modulo division is an 



y 



iterative algorithm: 
N „ R 



where 

N is the dividend, 
D is the divisor, 
Q is the quotient, 
20 R is the remainder, and 

0<R<D. 

Assuming Q is an n-bit positive integer, the value of Q can be denoted as a sum of 
binary positions: 

T>0*-| Q = ZQj-2 J (2) 

7=0 



* 2 5 where Q i is the jth bit of Q,Qj = {0, 1 } for all j 
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The traditional iterative algorithm begins by initializing R to N, i.e., R(n)=N and 
initializing Q n to 1 since N is a positive integer. Then the traditional algorithm calculates 
RO) and Q i for j=(n-l) to j=0 iteratively. The complete algorithm, for dividing an n-bit 
positive integer N by an n-bit positive integer D resulting in an n-bit positive quotient Q 
5 and n-bit positive remainder R, is shown in Figure 1 . 

An example of the traditional iterative algorithm is set forth below: 
Example 1: 



15 



10 



Assuming N=25, D=7, and n=3 => 25 mod 7 = 4 

stepAl: R(3)=25, Q 3 =l,j=2 

stepBl: j<0?No 

step CI: Q 3 =l?Yes 

step D 1 : R(2)=R(3)-D2 J =25-(7)(4)=-3 

step El: R(2)<0?Yes 

stepFl: Qj=Q 2 =0 

stepGl: j=j-l=2-l=l 

step B2: j<0 ? No 

step C2: Q 2 =l ? No 

step HI : R(l)=R(2)+D2 J =-3+(7)(2)=l 1 

stepE2: R(l)<0?No 



20 



step II: 
step G2: 
step B3: 
step C3: 
step D2: 
step E3: 
step 12: 
step G3: 
step B4: 
step Jl : 



Q.=i 

j=j-l=l-l=0 
j<0 ? No 
Q,=l ? Yes 

R(0)=R(1)-D2'=1 1-(7)(1)=4 
R(0)<0?No 

Qo=l 

j=j_l=0_i=_i 
j<0 ? Yes 

R = R(-1)=R(0)+(1-Q 0 )D=4+(1-1)7=4 



25 



Docket No. 2925-0272: 



According to the traditional iterative algorithm, as illustrated in Figure 1 and 
Example 1, the remainder of 25 mod 7 is correctly computed to be 4. 

The traditional iterative algorithm for determining the remainder of modulo 
division is a general purpose method, in that any two positive integers may be entered for 
5 N and D. However, as illustrated in Example 1, the traditional iterative algorithm requires 
numerous computations due to its iterative nature. For example, if D is an n-bit integer 
with a value D=2": 1 and n is any positive integer less than or equal to (D-l) 2 , i.e., n is a 2*n- 
bit integer at maximum, the traditional iterative method requires 6*n numbers of condition 
testing, 2 n numbers of multiplications (or shifts) and 2-n additions. Therefore, the total 
p 10 number of operations required is 10*n, excluding any value-assigning operations. 



technique for calculating the remainder of modulo division. The present invention requires 

1 5 significantly fewer operations than the traditional iterative technique for the same 

calculation. Furthermore, the number of calculations requires in the present invention is 
independent of the number of bits of the divisor in the modulo operation. 

Two requirements of the non-iterative technique of the present invention are that 
the value of the divisor D should be equal to 2 n -l (where n is the number of bits of the 

2 0 divisor D) and the value of the dividend N should be less than or equal to (D-l) 2 , but great 
than or equal to zero. 

Because the two constraints set forth above are basic constraints of Reed-Solomon 
coding, the present invention is extremely useful in applications that implement Reed- 
Solomon coding. Reed-Solomon coding involves algebraic operations in a Galois field. 

2 5 Reed-Solomon coding is a type of forward-error correcting coding that is used extensively 
in data communications. In Reed-Solomon coding, both conditions of the technique of the 
present invention are met and the algorithm of the present invention greatly improves the 
Reed-Solomon coding speed. In particular, the present invention is at least four times 
faster than the traditional iterative algorithm for a 16-bit fixed-point digital signal 

30 processor with special instructions supporting iterative division. Furthermore, if the 16-bit 



SUMMARY OF THE INVENTION 



The present invention is directed to a method and apparatus for calculating the 
remainder of a modulo division. The present invention is directed to a non-iterative 
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fixed-point DSP has no special division instruction, the algorithm of the present invention 
is at least twelve times faster than the traditional iterative algorithm. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 illustrates a flowchart of the traditional iterative method for determining 
the remainder of a modulo division. 

Figure 2 illustrates a hardware arrangement of the present invention of one 
embodiment. 

Figure 3 illustrates a flowchart of the technique of the present invention for 
determining the remainder of a modulo division. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

Figure 2 illustrates an apparatus of the present invention in one embodiment. In 
particular, Figure 2 illustrates a coder 10, which includes a processing unit 12, which 
implements the algorithm of the present invention. The processing unit 12, which performs 
the algorithm of the present invention, may receive a computer program to implement the 
algorithm of the present invention via an article of manufacture 14 or propagated signal 
20. The article of manufacture 14 further includes a medium 16, in addition to the 
computer program 18. 

The processing unit 12 may be any analog or digital processor, and either 
hardwired or software programmable to carry out the algorithm of the present invention. 
Processing unit 12 could be a general purpose processor, a digital processor (DSP), an 
algorithmic logic unit (ALU), or any other processing element, either discrete or 
integrated, which performs the algorithm to be described below. The coder 10 could be 
any type of coder which includes the processing unit 12, such as, for example, a Reed- 
Solomon coder. 

The processing unit 12 implements a non-iterative technique for calculating the 
remainder of modulo division. The number of calculations performed by the processing 
unit 12 is independent of the number of bits of the divisor in the modulo operation. Two 
requirements of the non-iterative algorithm of the present invention are that the value of 
the divisor D should be equal to 2 n -l (where n is the number of bits of the divisor D) and 
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the value of the dividend N should be less than or equal to (D-l) 2 , but greater than or equal 
to 0. If the two conditions set forth are met, the remainder R of M mod D is determined by 

"]~ b b summing the upper — and lower — bits of the dividend N to produce the remainder R. 

— 1 . — 1 

The two conditions described above enable the processing unit 12 to determine the 
5 remainder of a modulo division operation with significantly fewer operations than the 
traditional iterative technique for the same calculation. 

The algorithm of the present invention will now be described in more detail. 
For a given dividend N and divisor D, computing N/D gives a quotient Q and 
ggji remainder R. The relation between these quantities is: 

^ N R /q\ 



01 

n 



where 0<R<D 
Rearrange the equation above gives: 

N-QD+R 
Equation 4 can be written as 



(4) 



ill 



Q(D + l)+(R-Q) i/R>Q (5) 



T|jJ 0^2^ N \(0-l).(D + l) + [(*-e) + (£> + l)] else 



Equation 5 can be reduced to 

N=Q'-(D+1)+R' (6) 

where, 

0 '=l Q ifR - Q 

1(0-1) else (7) 



R ,J (R-Q) ifR^Q 



{[(R-Q) + (D + 1)] else 
2 0 Comparing Equation 4 and Equation 6 it can be seen that Q' and R' are the quotient and 
remainder of N(mod (D+l)) and Equation 7 gives the relationship between Q, R and Q\ 
R'. Adding Q' and R' gives: 



r^otr Q'+«'={ R R +D 



ifR^Q (8) 
else 
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Equation 8 shows the sum of Q' and R' can be equal to R or R + D. Given R is a 
nonnegative integer less than D as shown in Equation 3, R can be obtained by comparing 
the sum of R' and Q' with D: 



_ T D j Q'+R' ifQ'+R'<D 

Tptno R = \ Q ,+ R ,_„ else 



D 



5 If the quotient Q' and remainder R of N(mod(D+l)) are available, the remainder of N(mod 
D) can be obtained using Equation (9). 

Given the above conclusion, now the question is whether there is an efficient way to 
calculate the quotient and remainder of N(mod (D+l)). Assuming D=2 n -1, 

r 2«-l 
N = X! a i ' • 2 ' where ° Jis the J th bit of N,a J = {0,1} for all j 

7=0 1 

1 0 then (D+l) = 2 n . This fact results in a quick method of obtaining Q f and R\ Assuming 
N<(D-1) 2 means N is less than or equal to 2-n bits. The binary representation of N is: 



T EP CT7£_ N = CZa k+ „ ■ 2 k ) • 2" + ][>,. ■ V 

frt /=0 



(10) 



where 



1,672, S^ + „-2*<2''^ga y .2^<2" 



*V _ ^ ~» — (11) 



15 Comparing Equation 6 and Equation 10 and assuming 2 n =(D+l), the quotient Q' and 
remainder R in Equation 6 can be obtained as: 

a. -7* -where n ?c tho (vt _u lr\th hit nf AT n 

k+n 



k=0 



Q'=Ys a k+n ' 2 * wherea k+n isthe(n + k)thbitofN,a k+n = {0,1} for all k 



= • 2 7 where aj is the jth bit of N 9 a,. = {0,1} for all j 

Equation (12) means that the quotient Q' and remainder R of N(mod (D+l)) are the 
numbers made up by the high and low n bits of N, respectively. 

20 A flow chart of the present invention is illustrated in Figure 3. Example 2, set forth 

below, is explained in conjunction with the flowchart illustrated in Figure 3. 



7 
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Example 2: 

IfN=25andn=3, 
Condition 1: D=2 n -1=7 => True 
Tj 0&t° Condition 2: 0<N<(D-1) 2 



5 0<25<36 => True 

Assuming N=25, D=7, and n=3 => 25 mod 7=4. 

N=25 is represented as a 2n bit number (five bits with a leading zero added) 

Oil 1001 

1=3 step A': a = high bits of N -» a = 0 1 1 

p 10 stepB': b = low n bits of N ->• b = 001 

Jit . stepC: c = a + b= 100 

p step D': c < D ? Yes 



step E': N(mod D) = c = 1 00 = 4 

According to the technique of the present invention, the remainder of 25 mod 7 is 
if! 1 5 correctly computed to be 4. Comparing the present invention to the traditional iterative 
b; technique, the present invention only requires 5 operations for determining the remainder 

R, as illustrated in Example 2, whereas the traditional iterative technique, as illustrated in 
Example 1, requires 21 operations. 

Another example illustrates the additional processing performed in step F f of the 
2 0 present invention. 

Example 3: 

IfN=15andn=3, 

Condition 1 : D = 2 3 -l=7 => True 
Condition 2 : 0<N<(D- 1 ) 2 

0<15<36=>True 
Assuming N=15, D=7, and n=3 => 15 mod 7 = 1 

N=15 is represented as a 2n bit number (four bits with two leading zeroes added) 

001:111 



r,o 



Z- 
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step A 1 : a = high n bits of N -> a = 001 

stepB': b = lownbitsofN->b = 111 

stepC: c = a + b=1000 

fV°1 0 stepD': 1000<0111 ?No 

5 step F: c = c-D = 1000-01 1 1 = 0001 

step E': N (mod D) = c = 0001 = 1 

Again, according to the technique of the present invention, the remainder of 1 5 
mod 7 is correctly computed to be 1 . 

The invention being thus described, it will be obvious that the same may be varied 



!:-y 



10 in many ways. Such variations are not to be regarded as a departure from the spirit and 
scope of the invention, and all such modifications as would be obvious to one skilled in 
the art are intended to be included within the scope of the following claims. 



